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Abstract 

In field theories where a metastable false vacuum state arises 
as a result of radiative corrections, the calculation of the rate of 
false vacuum decay by bubble nucleation depends on the effective 
potential and the other functions that appear in the derivative 
expansion of the effective action. Beginning with the Nielsen 
identity, we derive a series of identities that govern the gauge 
dependence of these functions. Using these, we show, to leading 
nontrivial order, that even though these functions are individually 
gauge-dependent, one obtains a gauge-independent result for the 
bubble nucleation rate. Our formal arguments are complemented 
by explicit calculations for scalar electrodynamics in a class of 
gauges. 



1 Introduction 



In addition to the minimum energy true vacuum state, many quantum field theories have 
one or more metastable "false vacua" that can decay to the true vacuum by the nucleation 
of bubbles of the stable vacuum. Methods have been developed for calculating the rate of 
this process either at zero temperature [0 or at high temperature ||. However, these must 
be modified in the case of theories in which symmetry breaking arises as a result of radiative 
corrections 0. While a scheme for dealing with such cases (at zero temperature) has been 
developed [Q, it leads to an expression for the bubble nucleation rate that is not manifestly 
gauge-independent. In this paper we address the issue of this gauge dependence. 

The standard approach |T| to the calculation of the bubble nucleation rate at zero tem- 
perature is based on finding a "bounce" solution of the classical Euclidean field equations. 
The nucleation rate per unit volume T may be written as 



where B is the Euclidean action of the bounce solution and A is an expression involving 
functional determinants that is generally equal to a numerical factor of order unity times a 
dimensionful quantity determined by the characteristic mass scales of the theory. 

A problem arises if radiative corrections modify the vacuum structure of the theory. 
Theories in which this happens generally have no bounce solution; even if a bounce does 
exist, the nucleation rate calculation based on the bounce is not reliable. However, by 
integrating out certain fields at the outset, one can derive a modified algorithm [|J that can 
be applied to this situation. The results of this method are conveniently expressed in terms 
of the functions that appear in the derivative expansion 



of the Euclidean effective action. (The dots represent terms containing four or more deriva- 
tives; these do not enter the calculation at the order to which we work.) These functions 
can in turn be expanded in power series in the couplings. For example, in a gauge theory 
with weak (i.e., 0(e 4 )) scalar self-couplings, the effective potential is of order e 4 and may be 
written, using an obvious notation, as 



r = Ae 



-B 






V 



off 



(1.3) 



while 



Z = 1 + Z e 2 H . 



(1.4) 
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The first step in this approach is to use the leading approximation to the effective action 
to determine a bounce solution (f>b(x) through the equation 



□*=^-. (1.5) 

The desired nucleation rate is then given by 

T = A'e- {Bo+Bl) (1.6) 

where 



B = J d x 
turns out to be of order e -4 while 

B 1 = J d A x 



:i-7) 



Kf(^) + ^M^)(^) 2 



is of order e -2 . The calculation of the pre-exponential factor is much more complicated than 
in the standard case; in particular A' cannot be expressed solely in terms of the functions 
appearing in Eq. (|1.2|) . Nevertheless, one finds that, just as in the standard case, A' is equal 
to a numerical factor of order unity times a dimensionful factor determined by the mass 
scales of the theory. 

Like any physically measurable quantity, the nucleation rate should be gauge indepen- 
dent. Since the leading terms in the effective potential are gauge independent, there is no 
difficulty in this regard with respect to either B or the bounce solution itself. However, 
both of the functions that enter in B\ are known to depend on gauge. Our goal is to show 
that, nevertheless, these combine to give a gauge-independent contribution to the nucleation 
rate. Although we do not explicitly examine the prefactor A', we expect that our methods 
could be extended — albeit with considerably more technical complication — to show that 
it too is independent of gauge. 

Our approach is based on the Nielsen identity 0, which describes the gauge dependence 
of the effective action, and which has been used to show that gauge-independent physical 
quantities can be obtained from a gauge- dependent effective potential 0. In Sec. 2 we 
present a compact derivation of this identity, following the method of Kobes, Kunstatter, 
and Rebhan ||. However, the usual form of the identity is not quite sufficient for our 
purposes. Instead, what we need is a series of identities, each of which gives the gauge 
dependence of one of the functions appearing in the derivative expansion (|1.2j ). Although 
the identity for the effective potential is well-known, the remaining identities are, to our 
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knowledge, new. In Sec. 3 we derive these from the master identity and then use them to 
give a general proof of the gauge- independence of B\. To complement this formal proof, we 
have verified the relevant identities by explicit calculations for the case of scalar quantum 
electrodynamics in gauges. These calculations, which expand upon the work of Aitchison 
and Fraser |§, are described in Sec. 4. Section 5 contains some concluding comments. Some 
two-loop effective potential calculations are presented in an Appendix. 



2 The Nielsen identity 

In this section we use the method of Ref. |J to derive the Nielsen identity. We consider a 
gauge theory with fields denoted by <p{. The classical action 5* is invariant under a set of 
infinitesimal gauge transformations of the form 

5 g <fn = A% (2.1) 

where the Af are linear operators. ( We will henceforth suppress the index (3\ for scalar 
electrodynamics, which we examine in greatest detail, there is only a single gauge parameter 
6 in any case.) By choosing a gauge-fixing function F(<fti) and introducing Fadeev-Popov 
ghosts 1] and f/, we can write the generating functional of connected Green's functions as 



exp [iW(J, F)] = I [Vfc] [Drj] [Vfj] exp 

where 



U(F) + i UViy)^) 



(2.2) 



is invariant under the BRST transformations 

£brst <f>i = C Aj?7, <5 B rst fj = ~^~£ F ' ^brst ^7 = (2.4) 

with ( an arbitrary Grassman number. 

In presenting the derivation, it is convenient to adopt a compact notation where 



(0{4>)) = e~ m / [X?0MT^]0(0,7^)exp 



il + i (fyJ^Uv) 



(2.5) 



for any operator O. Now note that if O is linear in the ghost fields, its odd Grassman 
character leads to the vanishing of this quantity. In particular, 

{f]G} = (2.6) 
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for any functional G[(f)(x)]. Applying the BRST transformation ( |2.4| ) to this equation results 
in the identity 

(s BRST [fj(x) G(x)} + ifj(x)G(x) J d A yJ\y) SbrstUv)) = (2.7) 
which may be rewritten, using the anticommutivity of rj and fj, as 

- (±F(x) G{x) + ^)||A^)\ =-ij d 4 yj\y) (A iV (y) fj(x) G(x)) . (2.8) 

Now consider the effect of an infinitesimal change F — > F + AF in the gauge-fixing 
function. Recalling Eqs. ( p.2|) and ( |2.3|) , we see that the change in the generating functional 
W[J] is simply the integral over x of the left hand side of Eq. ( [2.8|) , with G set equal to AF. 
Hence, 

AW = -i I d A x d A y J\y) (A^(y) f)(x) AF(x)) . (2.9) 



Recalling that the effective action is related to W[J] by the Legendre transformation 

f a SW r A 

S cS = W - d 4 xJ\x) — — = W- d A xJ\x) <j>i(x) (2.10) 

we find that 

AS cS = z ( d 4 xd 4 y^^(A iV (y)f)(x)AF(x)) 1Pl (2.11) 

where the subscript 1PI indicates that only the contributions from one-particle irreducible 
graphs are to be included. In particular, an infinitesimal change d^ in the gauge parameter 
is equivalent to the choice AF = —(F/2£)d£. Hence, 

_dS e fi i f 4 4 5S e s 



^ = ~J d^xd'yj^iAMv) fj(x) F(x)) 1PI (2.12) 



d4 yjTT, H Am,y} (2-13) 
where 

Hj[<j>(z),y\ = - d 4 x{A i r 1 {y)n{x)F{x)) xm . (2.14) 
Eq. ( |2.13| ) is the Nielsen identity. 
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3 Derivative expansion of the Nielsen identity and for- 
mal proof of gauge independence of the nucleation 
rate 

To study the gauge-dependence of the bubble nucleation rate, we will need a set of identities 
that are obtained by making derivative expansions of both sides of Eq. ( 12.131) . For simplicity, 
consider the case where the effective action depends on only a single field 4>(x). There is 
then only a single functional H[4>(x),y], which can be expanded as 



H[^x),y]=C{<f>) + D^){d ll <f>) z + 



(3.i; 



where all terms on the right are understood to be evaluated at point y and the dots represent 
terms with more than two derivatives. Inserting this, together with the expansion (|1.2|) of 
the effective action, into Eq. ( |2.13| ), gives 



d_ 



d x 



V cff (0) + -Z(0)(^0) 2 + 



d x 



C(cf)) + D(<f>)(d,<f>) 2 + 



dV 



cff 



1 r)7 

+ - d, [Z{4>) d,<j>\ + ■ ■ 



(3.2) 



If this identity is to hold for arbitrary <f)(x), then not only must the integrands on the 
two sides be equal point by point, but the terms with equal number of derivatives must be 
separately equal. Thus, the terms with no derivatives obey 



dV 



cff 



C 



dV 



cff 



while from the terms with two derivatives we obtain 

C dZ -C 9Z I 2D dVCS I 2Z d ° 

d£ d(j) dcj) dcj) 



(3.3) 



(3.4) 



(Eq. ( |3.3|) , which can be obtained immediately from Eq. (|2.13|) by choosing <p(x) to be a 
constant, appears in Ref. 

We now specialize to the case of a gauge theory with gauge coupling e and scalar self- 
couplings of order e 4 . As indicated in Eq. (|1.3|), the effective potential begins with terms of 
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order e 4 , while Z((j>) = 1 + 0(e 2 ). Analysis of the relevant graphs shows that C(0) starts at 
order e 2 and D((j>) is of order unity. The terms of order e 4 in Eq. ( |3.3|) yield 



dVf 



0. 



(3.5) 



Now recall that the bounce solution <f>b(x) is determined, through Eq. (|1.5|) , by V^F. Since 
Eq. (|3~5|) shows that the latter is gauge independent, both <f>b(x) and B , the leading contri- 
bution to the exponent of the nucleation, are independent of £. 

To study the gauge dependence of B x , we need the order e 6 terms of Eq. (|3.3|) , 



dVf 



a 



dVf 



(3.6) 



as well as the terms of order e 2 in Eq. 

1 dZ e 2 _ 8C e 2 

These equations, together with Eq. ( [L.8|) , imply that 

dB x 



(3.7) 



4/ 


d 4 x 


/a 




J d 4 x 


C e 2 





Vf + -Z e2 (dM 2 



dVf dC e 

-7T— + 



cff 



+ {d tl C e 2){d^ h ) 



cff 



□t 



(3.8) 



where all quantities are to be evaluated with <p(x) set equal to the bounce solution <fib( x )- 
Eq. (|1.5|) , which determines the bounce, shows that the last expression on the right hand 
side must vanish, and hence that 



(3.9) 



This verifies that, at least up to pre-exponential terms of order unity, the bubble nucleation 
rate is gauge independent. 
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4 Scalar Electrodynamics 



4.1 Basics 

We now illustrate these formal arguments by explicit calculations for the case of scalar 
electrodynamics. The Lagrangian, which we write in the form 

£ = + \( d & ~ e/ W 2 + ~(^$2 + eA^f - vm (4.1) 

with 

V(<f>) = im 2 $ 2 + (4.2) 

and $ = ($ 2 + c^ 2 ,) 1 / 2 , is invariant under the gauge transformation 

6 g A^ = d^9, Sg®! = e^ 2 6, 5 g $ 2 = -eQtf . (4.3) 

If m 2 > 0, the tree-level potential has a minimum at $ = 0. In order that one-loop effects 
be able to change the vacuum structure and give a symmetry-breaking minimum at $ = 
($) 7^ 0, we must require that both A and m 2 be anomalously small, of order e 4 and e 2 ($) 2 , 
respectively. 

For calculating the bubble nucleation rate it is sufficient to evaluate the terms in the 
derivative expansion of the effective action for <3> 2 — A^ = 0. With this in mind, we will 
consider the class of gauges determined by the gauge-fixing function 

F = (d^ + ev^ 2 ) . (4.4) 

(The gauge-dependence of the effective potential in these gauges was studied in detail by 
Aitchison and Fraser ||; in the following discussion we will make use of a number of their 
results.) The Nielsen identity (|2.13|) then involves only the single functional 

H 9l [<f>(z),y] = -- J d 4 x(($ 2 (y) rj(y) fj(x) (d^(x) + ev<S> 2 (x))) . (4.5) 

The effective action can be obtained as the sum of one-particle irreducible vacuum graphs 
in the theory obtained from the Lagrangian ( |4.1| ) by making the shift $i — ► $i + <fi and then 
dropping all terms linear in the quantum fields. The vertex factors for these graphs can be 
simply read off from the resulting Lagrangian in the standard fashion (see, e.g., Ref. ||). 
The propagators require a bit more work. Following the usual approach, one would obtain 
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from the Lagrangian (together with the gauge-fixing and ghost terms) the effective $1, $2, 
Ap, and ghost propagators 

G '<*> = w^kw (4 - 6) 



G^(fc) = G^(fc) + G^(fc) 



-9,u + ^ i[((P-^)-eV]^ 



A; 2 - e 2 2 D(Jfe) 



(4i 



G 3 = 12 ■ 2 A ( 4 - 9 ) 



as well as the mixed $2 -An propagator 



G^(k) = cw ;: ft " (4.10) 



e( £0 + t>) kf, 
D(k) 

where the momentum flow is understood to flow from the $2 end to the A^ end. In these 
expressions 

D(k) = k 4 - k 2 (ml - 2e 2 v<t>) + e 2 (j) 2 (e 2 v 2 + £m 2 ) (4.11) 

and 

m 2 (0)= m 2 + ^ 2 = y"(0) (4.12) 

^(0)=m 2 + V = ^. (4.13) 

These propagators are not quite what we need. Our assumption that A is of order e 4 
not only makes some of the one-loop terms comparable to the tree-level terms, but also 
implies that some multi-loop graphs are not suppressed relative to graphs with fewer loops; 
specifically, the insertion of transverse photon loops along a scalar propagator does not 
increase the order of the graph. To restore the validity of our expansion, these insertions 
must be summed. This can be done simply by replacing the propagators given in Eqs. (fO|- 
by "dressed" propagators in which the fh 2 are replaced by 

m 2 (0) = V^<j>) (4.14) 
m 2 M = W l - (4-15) 
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(To avoid double-counting, subtractions are needed for certain graphs with two or more loops; 
these corrections only affect contributions of higher order than those we will be considering. ) 

Before we proceed to verify the identities, there is one more issue to be addressed. Many 
of the graphs contributing to the effective action have divergences that must be cancelled 
by appropriate counterterms. We will not display these explicitly, but all divergent integrals 
should be understood to be made finite by some gauge-invariant renormalization scheme 
(e.g., minimal subtraction in the context of dimensional regularization); when we refer to 
the magnitude of an integral, this should be understood as referring to the magnitude of its 
finite part. 



4.2 The identity for the effective potential 

To order e 2 the function C(4>) entering the identity ( |3.3| ) receives contributions only from 
the two graphs shown in Fig. 1. These combine to give [] 



C P 



ie r d A k 



2 J (2tt) 4 (k 2 + e 2 vcj)) D(k) 
ie 2 (j)£ r d 4 k 1 



e(£0 + v)k- ev(k 2 - f ) 



(2tt) 4 £>(£;) 



(4.16) 



The effective potential is obtained by summing the graphs with vanishing external mo- 
menta. The one- loop contributions may be split into three parts. First, the graphs with a 
transverse photon loop give a contribution 

3/ '' <rk ln(k 2 -e 2 <p 2 ). (4.17) 



2 J (2tt) 4 

From dimensional arguments, this integral is clearly of order e 4 4 . Hence, it combines with 
the tree-level potential to give 

Vf = \m 2 <P 2 + ^ 4 - | / Hk 2 ~ eV) • (4.18) 



This is manifestly gauge-independent, in accordance with Eq. (p.5|). A second gauge- 
independent contribution, coming from the graphs with $i-loop graphs, is 

' ' ,/,/,: ln(7,- 2 -/,/;'). (4.19) 



2 J (2vr) 4 



1 Apart from an overall sign arising from a difference in the definition in C(</>), these 
expressions are the same as those appearing in Ref. B 
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This is of order mf, and hence contributes only to r. Finally, there is an order e contri- 
bution 

^efr i f d k 



J \hiD(k) - 2 ln(fc 2 + e 2 ^)] (4.20) 



l-loop 2 J Q^y 

in which the first term arises from graphs with $ 2) longitudinal photon, or mixed scalar- 
photon propagators while the second is due to those with a single ghost loop. 

In addition to these one-loop contributions, there are a number of two-loop graphs that 
contribute to . Some two-loop graphs have already been included in Eq. (|4.20| ) as a 
result of the replacement fhf — > m 2 , including in particular the ^-dependent "figure-eight" 
graph with one transverse photon loop and one $2 loop. In the Appendix we show that 
although a number of the remaining graphs are ^-dependent, they add togetheif] to give a 
gauge-independent contribution to . Hence, we can combine Eq. ( |4.20| ) with Eqs. ( f4.11|) , 
(gap , and ( ETI6D , to obtain 

e ^T = e =Ce2 ^> (4 ' 21) 



thus verifying Eq. (|3.6|). 

4.3 The identity for Z(<f>) 



We now turn to the identity (|3.4j ), which we will verify to order e 2 . We begin by recalling that 
Z((j)) can be calculated from the sum of one-particle irreducible graphs with one external line 
carrying momentum p, another carrying momentum — p, and all others with zero momentum. 
If the contribution of graph j is denoted by Ij{p 2 ), then 



-yi 



(4.22) 

p 2 =0 



Although there are many one-loop graphs contributing to the scalar self-energy, we will 
need to calculate only a few. Those graphs with quartic vertices are independent of the 
external momentum and hence do not contribute to Z((p). Because of our assumption that 

2 The cancellation of the gauge-dependence among these graphs can be understood 
by considering the case A ~ e 2 , where the loop expansion is completely equivalent to 
an expansion in e 2 . Apart from the appearance of m 2 rather than m|, the one-loop 
approximation to C(4>) is precisely the same as C e 2 of Eq. ( 4.1 6| ) . Hence, Eq. fl3.3| ) can be 



satisfied both at the one- loop level in that case (as was shown in ||) and in at 0(e 6 ) in 
our case only if this cancellation among the two-loop graphs occurs. 
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A = 0(e 4 ), all graphs with a vertex arising from the scalar self- interaction are at least of 
order e 4 and can also be ignored here. Finally, the self-energy graph with a single ghost loop, 
although of order e 2 , is {-independent. Thus, the entire {-dependence of Z e 2(4>) comes from 
the four graphs shown in Fig. 2. 

It is convenient to consider separately the terms containing the transverse and the longi- 
tudinal parts of the photon propagators. The only contribution with two transverse propaga- 
tors, from graph b, is manifestly independent of { and so can be neglected here. Graphs a, b, 
and c each give contributions with a single transverse propagator; although separately these 
each contain {-dependent terms of order e 2 , their sum is easily seen to be {-independent to 
this order. 

This leaves the terms containing only longitudinal photons. These may be written as 
where 

b a = _[(p + kf - {e 2 2 ] ({A; 2 - eV) [k ■ (2p + k)} 2 k~ 2 

b b = 2e 2 2 ({A; 2 - eV) [{(p + kf - eV] [k ■ (p + A;)] 2 A;' 2 (p + A;)" 2 

b c = -4e 2 2 ({0 + v) ({A; 2 - e 2 v 2 )[k ■ (p + k)} [k ■ (2p + k)} k~ 2 

b d = e 2 ({0 + v) 2 (k 2 -p 2 )k-(k + 2p) . (4.24) 

(In these expressions we have omitted terms proportional to m\ since the contribution of 
these is at least 0(e 4 ). Note also that we have included a factor of two in the the contribution 
from graph c to take into account the fact that reversing the direction of p gives a second 
graph with the same value.) Summing these expressions, we obtain 



b a + b b + b c + b d = {(p 2 - 2p • k - k 2 )D{p + k) - 2{p 2 (e 2 (f)v + e 4 : 



V 



+ 



+ • • • (4.25) 

where the dots denote terms that are either {-independent, of order p 3 , or else proportional 
to m 2 , and thus of higher order in e. Because of their antisymmetry under the interchange 
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k 2 < — > (p + k) 2 , the contributions of the terms in curly brackets to the integral in Eq. ( |4.23| ) 
cancel. Inserting the remaining terms into the integral and keeping terms proportional to 
p 2 , we find that 



Z e 2 




d 4 k 1 2 (e 2 (f)v + e 4 (j) 2 v 2 ) 



+ £— independent terms 



(2tt) 4 [D(k) [D(k)} 2 
d 4 k [ 1 dD(k) 



+ £— independent terms 



(2tt) 4 [£(*;) [£(£:)] 2 <90 



90 



+ £— independent terms . 



(4.26) 



(In going from the first to the second line, contributions proportional to drn^/dcj) have been 
neglected as being of higher order.) 

Differentiating this with respect to the gauge-parameter £ gives 



This verifies the identity ([3.71) . 

5 Concluding remarks 

In this paper we have have shown how the Nielsen identity that describes the gauge depen- 
dence of the effective action can be converted into an infinite series of identities, one for each 
of the functions appearing in the derivative expansion of the effective action. Using these 
identities, we have shown, to leading nontrivial order, that one obtains a gauge-independent 
result for the bubble nucleation rate even in theories where the calculation of this rate in- 
volves the gauge-dependent higher order contributions to the effective action. This provides 
one more example to show that the gauge- dependence of the effective action does not prevent 
it from being a useful tool for obtaining gauge-independent physical results. 

As an explicit example, we have verified the identities for the ^-dependence of V eS (<fi) 
and Z(cj)) in the class of gauges defined by the gauge-fixing function ( |4.4j ). In fact, these 
gauges actually depend on a second parameter, v. (Note that nothing in our calculations 
requires that v be equal to the vacuum expectation value of 0.) Working from Eq. ( |2.11| ), 





<90 



+ 0(e 4 ). 



(4.27) 
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we find that 

dV eS dV eS 

(5 - 1) 

and 

Mc.oz +2I r?£L + 2Z «?L. (5.2) 

dv dtp o<fi o<p 

where C v (4>) and D v ((j)) are obtained from the derivative expansion of 

Hl^{z),y] = -ie 2 v j d A x((<$> 2 (y) r,(y) fj(x) $ 2 {x)) . (5.3) 

(Note the absence of the factor of 1/2 relative to Eq. ( |4.5| ).) In particular, the leading 
contributions to these identities comes from 

. 2 r d*k (k 2 - £e 2 2 ) K A , 

^ = - ieV ) (2^ (P + eV)DW (5 - 4) 

The Fermi gauges, defined by F = <9 M A M , can be obtained from the gauges we have 
considered by setting v — 0. However, because of the infrared divergences that afflict these 
gauges, the limit v — > is somewhat nontrivial and the verification of the identities for these 
gauges must be done separately |fTC|| . To see the cause of these difficulties, note that our 
assumptions about the magnitude of A and m 2 imply that if v = the zeros of D(k) occur at 
values of k 2 of order e 4 2 . This has the effect of making some two- loop graphs (beyond those 
resummed by the conversion of m 2 to m 2 ) comparable to one-loop graphs. For example, in 
the calculation of the quantity dC e 2/d(p on the left hand side of the Nielsen identity for Z e 2, 
the terms involving dm 2 /d<j) are no longer higher order. The corresponding terms on the 
right hand side of the identity come from contributions to Z e 2 due to two-loop graphs and 
one-loop graphs with vertices proportional to A; both types of contributions can be neglected 
for generic nonzero values of v. 



A Appendix 

In this appendix we show that, although individual two- loop graphs give gauge- dependent 
contributions to V^f, their sum is ^-independent. The first step is to identify the relevant 
graphs. All graphs with vertices proportional to A give higher order contributions, and so can 
be omitted. Similarly, any graph with a loop containing only $i propagators is proportional 
to a power of m\ and hence of higher order. Finally, there is a two-loop graph containing a 
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ghost loop, but it is manifestly gauge-independent. The only nonzero graphs remaining are 
shown in Fig. 3. 

It is convenient to decompose the photon propagators into transverse and longitudinal 
parts, and to examine separately the contributions from each. Consider first the contributions 
involving transverse photon propagators. The part of graph a involving two such propagators 
and the part of graph e involving one transverse photon have already been included in the 
one- loop calculation by the resummation that converted the m 2 to the m 2 , and hence should 
be omitted. This leaves the portions of graphs a, b, and d that involve only a single transverse 
photon each. The contribution of these to the effective potential is 

_ 2 f _d^^p 9^-^ 

6 J (2tt) 4 (2tt) 4 (Jfe 2 - e 2 2 ) D{p) [(k + pf 



mf] 



2„.2\ 



2eV (£p 2 - e 2 v 



(k + 2p) fl (k + 2p) u - 2e 2 0(£0 + v)(k + 2p) il p v 



= 2e 2 



d A k d 4 p 



1 



mf] 



1 - 



(P • kf 
k 2 p 2 



(A.l) 



(2tt) 4 (2tt) 4 (k 2 -e 2 2 ) [(k + pf 

where the dots represent terms proportional to m\. Not only is this result independent of £, 
but examination of the integrals shows it to be in fact of order e 8 . 

This leaves us with the terms involving only longitudinal photons. Let us denote the 
corresponding contribution from graph j by Jj. For the first four graphs this may be written 
in the form 

d*k d A p 1 



aj(k,p) 



(2tt) 4 (2tt) 4 [(jfe + pf - m 2 ] D{k) D{p) 
Omitting terms proportional to m 2 , whose effects are of higher order, one finds that 

a a = e 2 2 (£k 2 - eV) {£p 2 - eV) (k ■ pf AT 2 p~ 2 
a b = -\ {ik 2 - eV) (p 2 - ^e 2 2 ) (k 2 + 2k • pf AT 2 

a c = ^ + vf (k 2 + 2k-p) (p 2 + 2k-p) 



(A.2) 



a d = -2e 2 (^0 + v) {£k 2 - e 2 v 2 ) (k +2k-p){k-p)k 



(A.3) 



The sum of these is 



a a + a b + a c + a d 



--(p + kf [(k 2 



(ip 2 - eV) + (k ■ p) [2p 2 i - e 2 <p 2 e" 



+| D{p) [ip + A;) 2 - p 2 ] + A(k, p) + 



(A.4) 
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where A(k,p) is an antisymmetric function of k and p and the dots represent terms that are 
either proportional to m|, and thus of higher order, or else ^-independent. 

When this sum is inserted back into the integral, the term containing A(k,p) vanishes 
because of its antisymmetry. The remaining terms give 



J a + Jb + Jc + Jd 



e 
~2 



I 



d 4 k d 4 p 



(2tt) 4 (2tt) z 



1 + 



[(k + p) 2 — m 2 ] 



(k 2 - ee 2 2 ) - eV) + (k.p) (2p\ - e 2 2 ^ 
e 2 t\ r d 4 k d 4 p [{p + k) 2 -p 2 } 



D(k)D(p) 



r dV 
J CM 



2 J (2tt) 4 (2tt) 4 [(k + p) 2 -m 2 ]D(k) 



+ 



(A.5) 



The terms in the first integral that are proportional to m\ are at least 0(e 8 ) and can be 
omitted. In the second integral, let us make the change of variable p — > p — k. The resulting 
integral is then clearly the product of two integrals, one of which is proportional to m 2 , and 
is thus also higher order. Hence, 



J a + Jb + Jc + Jd 



e 
~2 



d 4 k {k 2 -ie 



(2tt) 4 D(k) 



d p (£p - e 2 v ) 



(2tt) 4 D(p) 



+ 



(A.6) 



This last expression is precisely equal to — J e . Hence, the two-loop contribution to VjF is 
^-independent, as was claimed. 
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Figure 1: The two graphs that contribute to C e 2. Photon, <£> 2 , and ghost propagators are 
indicated by wiggly, long-dashed, and short-dashed lines, respectively. 
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